Point-like spin-dependent interaction 
in calculations of self-energy ladder diagrams 
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An instantaneous and zero-range spin-dependent interaction, derived by summing an infinite 
number of electron-fiole ladder diagrams witliin a local approximation, is analyzed as a function 
of the electron gas density and the relative spin polarization. The strength of such an interaction 
is defined by an integral of a statically screened Coulomb interaction with a spatially localized 
weighting factor. This weighting factor represents the mutually uncorrelated motion of an electron- 
hole pair in singlet or triplet spin states (5*^ = 0, ±1). An implementation, based on a Yukawa- type 
interaction with a spin-polarization-dependent Thomas-Fermi screening length, is given. 
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A rigorous study, based on quantum many-body per- 
turbation theorjii as applied to excited quasiparticlc 
states in itinerant ferromagnetic materials, requires an 
approach to the quasiparticlc self-energy accounting 
for the contribution of both charge- and spin-density 
fluctuationsi^ Within this approach, due to the corre- 
spondence between multiple electron-hole (e — h) scat- 
tering events and a spin fluctuation)2i4 the self-energy 
can be treated as an integral over the four-point e — h 
scattering amplitudei^ 

The theoretical problem, which appears already at the 
homogeneous electron gas level, lies in this amplitude and 
can be largely simplified in the following ways: (i) by us- 
ing a point-like interaction with strength defined by the 
Hubbard U parameterji (ii) by introducing a local pseu- 
dopotential obtained within some self-consistently closed 
procedurejii (iii) by reducing an infinite set of diagrams 
to a local-field factoid tabulated and parameterized by 
making use of accurate quantum Monte Carlo data for 
the electron gas. 

Apart from the obvious theoretical interest for a pro- 
totype homogeneous system, real materials give an ad- 
ditional important challenge. In evaluating the quanti- 
ties listed above, one should properly take into consid- 
eration the band-structure effects which surely influence 
calculated characteristics of excited quasiparticlc states. 
However, the obtained results are very sensitive to an 
approximated form chosen for these quantities There- 
fore, the question of how in ab initio calculations to sum 
an infinite number of diagrams, which contribute to the 
four-point e — h scattering amplitude, is still open. 

The standard formulationiS for excited states starts by 
defining a dynamical, spin-independent dielectric screen- 
ing of the bare Coulomb interaction between two test 
charges. This concept of a polarizable medium allows 
a straightforward and successful implementation of the 
Hartree-Fock-like technique with such a screened inter- 



action instead of the bare one (the GW method), since 
it is based on a complete set of single-particle states at 
the Hartree (mean-field) level. Thus the required ap- 
proach naturally should start with the above screened 
interaction to consider multiple scattering effects in the 
important e — h channel. 

Here we present a theoretical study which rests on 
Hedin's conccptiiS We derive a point-like spin- dependent 
interaction. The foremost advantage of using the latter 
is that a number of many-body problems becomes more 
easily solvable. To achieve this goal we shall use the 
approach of Ref. ^3 which was formulated for uniform 
systems in momentum space. This approach is based on 
a treatment of the self-energy ladder diagrams and does 
not contain adjustable parameters or quantities defined 
outside the scope of the approach. A key quantity of 
the approach is a local interaction Wa-a' that (i) ensures 
a simplification similar to Hubbard models, (ii) includes 
charge and spin fluctuations, and (iii) has a transpar- 
ent connection with the exchange part of the many-body 
local field factor. 

Accomplishing our goal, first, we perform a real-space 
analysis of the approach of Ref. [Ill Such an analysis 
allows one to derive basic equations of the approach to 
be applied to a broader spectrum of materials. Unless 
stated otherwise, atomic units, = fi, = m = 1, arc used 
throughout this work. By applying standard notation for 
the coordinate, 1 = (cri, ri, ti), in which the quasiparticlc 
spin is denoted as cri , the Bethe-Salpeter equation defin- 
ing the e — h scattering amplitude (the T matrix) can be 
written as^ 



T(l,2|3,4) = T°(l,2|l',2') [(5(l'-3)(5(2'-4) 

+ A'°(1',2"'|3',4')T(3',4'|3,4)] , (1) 

where T°(l, 2|1', 2') is the irreducible e — h vertex and the 
free propagator for e — h pair is given by the product of 
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the Green functions as K°{1,2\1' ,2') = iG{l,l')G{2' ,2). 
In Eq. , an integral over space-time coordinates and a 
sum over spin are implied for repeated variables. 

To perform a selective summation of ladder diagramsii 
along the line of Hedin's expansiorkiS for the self-energy 
in terms of a dynamically screened Coulomb interaction 
W, the irreducible e — h vertex can be chosen as 

rO(l, 2|3, 4) = Wil, 2)S{1 - 3)<5(2 - A)S^,.J,,.„ 

where in the right-hand side the spin variables are explic- 
itly written out and space-time coordinates are abbrevi- 
ated as 1 = (ri,ti). As a result, the T matrix is defined 
by the ladder approximation to the Bethe-Salpeter equa- 
tion as 

T,,,(l,2|3,4) = W{l,2)[5il-3)Si2-A) 
+ J dl'd2'X0,,(l,2|l',2') 

X T,,,(l',2'|3,4)], (2) 

where the e ^ h propagator expressed in terms of the 
Green functions diagonal in spin space. In Eq. (O, there 
is no sum over spin implied and all integrals are written 
explicitly. The T matrix (|2J) describes propagation of a 
mutually correlated e~h pair carrying spin = 0, ±1^ 
This correlation is realized by repeatedly interacting pair- 
constituents via 1^(1, 2) in the intermediate state. 

In real space, the approximation done in Ref. ITll results 
in the T matrix given by 

T,,.(l, 2|3,4) = f,,,(l|4),5(l - 2)6(3 - 4), (3) 

where 

f,,.(l|4) = W^a'ilA) + J dl'd2'W^,,{l,r) 

X X°,,(lM'|2',2')f,,,(2'|4). (4) 
The interaction Waa'i^, 4) is defined by the equation 

rfl'rf2'if°,,(l,l|l',l')W?..'(l',2')i^°.'(2',2'|4,4) 

I dl'd2'X0,,(l,l|l',2')Ty(l',2')ifL'(l',2'|4,4). (5) 

The obtained T matrix is a local one describing scattering 
processes in which the coordinates of the e — h pair both 
for initial states and for final states coincide. This local- 
ity essentially simplifies the self-energy evaluation as an 
integral over the T matrix. Actually, owing to Eqs. (O 
and Q), the T-matrix contribution as an additional term 
to the GW self-energy 2) = iG^{l, 2)W{1, 2) can 

be cast into the GVF-like form given by 



E^(l,2) = -z^G,,(l,2)r,,,(l|2), 



(6) 



where 



r.^,(l|2) = / dl'd2'iy<,,,(l,l')A'L'(lM'|2',2') 



f,,,(2'|2)-iy,,,(2',2) 



(7) 



Thus, in real space we have obtained the T-matrix contri- 
bution by exploiting the local approximation which seems 
to be reasonable for the corrections to the GW term (see, 
e.g., Ref.m. 

Now we can address to our main goal, deriving an in- 
stantaneous and zero-range potential 



iy,,Ki,2)-v;V<5(i-2). 



(8) 



It is obvious that in this case the irreducible e — h vertex 
becomes highly local as in Hubbard models. Further, we 
will be guided by the results of Ref. where for a uni- 
form system and at small four-momentum transfer along 
the e — h channel the local interaction W^a' was obtained 
to be constant in momentum space. Consequently, its 
real space equivalent corresponds to a 5-function interac- 
tion we need: the resulting T matrix has the form which 
can be obtained from Eq. (|2Jl by using the zero-range 
spin- dependent interaction (jS]) instead of M^(l,2). 
In terms of the noninteracting e — h Green function 



[I - (eka')] nF(ek-i-q^) 

- Ck+qcr + Cko-' - "7 

'•^f (eka-') [1 - »-F(ek+qg)] 

^ ~ Ck+qcr + CkfT' + ill 



where np is the Fermi distribution function and the en- 
ergy eko- is measured from the Fermi energy ep, the free 
e — h propagator is given by 



(2^ 



We introduce a quantity which characterize the mutually 
uncorrelatcd motion of the e — h pair as 



1 



iC°,,(0,0)7 (2^) 



limgL.'(q,0). (9) 



q^O 



This quantity is unity at r = and tends to zero as ^ 1/r 
at r — > oo. Thus, we can write the point-like interaction 
strength in real space as 



dv\W{r)?V{v), 



(10) 



with a statically screened Coulomb interaction denoted 
as V{t) = W{v,bj = 0). In this expression, [^^.^./(r)^ can 
be considered as a spatially localized weighting factor. 

In models based on the underlying Hartree approach, 
the interaction in Eq. (|10|1 can be approximated by 
a simple Yukawa form (1/r) exp(— Ar), where A speci- 
fies a characteristic length of the screening in the sys- 
tem [e.g., in Thomas- Fermi approximation A = Xtf = 
[Akp /t^Y^"^]. For completeness, we have to note that in 
Hartree-Fock-type approximations for electron-electron 
interactions the screening length cani^iiStiii^ differ (A ~ 
kp) from the Thomas-Fermi one. Hedin's expansion, be- 
ing our basic frame in the present work, rests on VF(1, 2) 
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FIG. 1: The spin-diagonal part of the point-like interaction 
strength V^cr over its value in the paramagnetic state as a 
function of at the SP C = 0.33, 0.67, and 1.0. 



which corresponds to the dielectric (Hartree) screening 
between two test charges. 

In order to highlight the physical meaning of the quan- 
tities C<T<T'(r) and V^^, entering Eq. |[TU|l . first, we con- 
sider the homogenous electron gas (HEG) in the para- 
magnetic state. In this case ^cro-'(r) [denoted as ^p(r)] 
can be easily found as Cp(r) = jo{kpr). Here ji is 
the spherical Bessel function of order / and kp is the 
Fermi wave vector. Thus, the point-like interaction of 
strength can be identified with the Fermi-Huang 
pseudopotentialiS — 27r(tan(5o//ifc)(5(r) at fc = fc^ is 
the reduced mass of the interacting pair) with the phase 
shift do evaluated within the first-order Born approxima- 
tion as tan^o ~ — (/iA:/27r) J dr[jo{kr)]'^V{r). 

As a digression, we would like to notice here that a 
somewhat similar idea was appliedSS earlier to the Lan- 
dau theory of liquid ^He. In Ref. [20l the quasiparti- 
cle scattering amplitude was described by the Fourier 
transform J e'^'^'^Vef /{r, k)dr of the effective potential 
Veff{r,k) = Va{r) + (l/3)fcVT/2(r), where Vq and V2 
are some local potentials. This effective-potential form 
addressed the question of non-locality at a given scatter- 
ing momentum k. In the forward direction (q = 0) and 
k = kp, this amplitude is formally very close to the in- 
teraction strength in question, because for small kpr one 
gets Mkpr)Y « 1 - (l/3)fc|,r2. 

Using the Yukawa form for V{y), one obtains the fol- 
lowing expression for the strength 



K = -^111(1 + 44^ 



TF 



(11) 



that in a week scattering regime reproduces the well- 
known expression for the interaction energy^ of the quasi- 
particles both of which are on the Fermi surface. This 
is inversely proportional to showing an expected be- 
havior since the scattering has vanishing effect by grow- 
ing density. 

For a spin-polarized HEG, when the electron and the 
hole belong to the same a subsystem [S^ = 0), ^^^(1") = 
jo{kpr), where the Fermi wave vector kp for spin a is ex- 
pressed asSi kp = kp{l — fT^y/^ with the paramagnetic 
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FIG. 2: The spin-non-diagonal part of ^o-ct' as a function of r 
at different values of the SP C,. 
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FIG. 3: The spin-non-diagonal part of the point-like interac- 
tion strength V^^ over its value in the paramagnetic state as 
a function of Ts at the SP C, — 0.33, 0.67, and 1.0. Inset: the 
ratio Atf(C)/-^tf as a function of (. 



value of kp. Here <t = +1 for spin- up (|) and cr = — 1 
for spin-down (|), respectively. The relative spin polar- 
ization (SP) C is given by ^ = |7i| — nj/n, where ricr is 
the spin a electron density, n = + being the total 
electron density. 

In addition to (r) , the C dependence of the strength 

V^^ comes from Atf(C) 



Atf^5E.(1-^C)i/3 Ob- 
tained at the long-wave limit from the standard ex- 
pression for the irreducible polarizabilitjiSiS^ in the spin- 
polarized case. As a net result, V^^ has the same form 
as Eq. Hll|l but with kp and Xtf{C) instead of fc^ and 
Xtf, respectively. 

Fig. n illustrates the spin-diagonal part of the inter- 
action strength plotted as the ratio V^^/Vp for rang- 
ing from 2 to 6 (the metallic density range) at differ- 
ent C parameters. Due to the nontrivial dependence 
resulting from £,ii{r) and ^(r) with Atf(C) in the ba- 
sic Eq. as the SP increases Vjj becomes greater or 
smaller than its value in the paramagnetic state at fixed 
Vg. On the other hand, the results for V^^/Vp show that 
at fixed the ratio demonstrates a monotonic increase 
with increasing This increase becomes smaller as 
increases. However, at C ^ 1 limit this monotonic behav- 
ior together with decreasing n| can lead to a nontrivial 
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C dependence of the corresponding contribution to the 
quasiparticle self-energy. 

The spin-non-diagonal part of ^aa' (r) corresponds to 
the case when the electron and the hole belong to differ- 
ent spin subsystems {S^ = ±1). This part is given by 

^TiW = (l - ^) bo(fcF^) + J2(fcjr)]. (12) 

Note that ^Ti(r) = ?iT(r) and at C = 0, Cnlr) = Cp{r). 
At the C ^ 1 limit, ^-rj,(r) tends to jo{kpr) + j2{kpr). 

Fig. 121 gives insight into the detailed (radial) behavior 
of the spin-non-diagonal part of the weighting factor ap- 
pearing in Eq. (|10|l . Multiplied by [r/r^]^, this factor has 
the first maximum at r O.Sr^ (O.SGr^ for C = 1) and 
the first node at r ~ l.Trg (LSGr^ for C — 1). Practi- 
cally, owing to the r-depcndence of the applied Yukawa 
interaction, one can focus on the part of ^ti('^) which 
corresponds to the r-range from zero to the mentioned 
node. Within this interval the considered function shows 
the weak (^-dependence. As a result, the behavior of 
plotted in Fig. |3| is mainly caused by the ^-dependence 
of the screening length Atf(C) (see the inset of the fig- 
ure). As is evident from the figure, for ( up to 0.67 the 
shown ratio V-^^^/Vp changes within ~ 5%. This means 
that Hubbard's idea on constancy is closely satisfied. 

A comparison of Figs. ^ and O unambiguously shows 



that it is V|| in the spin-diagonal part of the point- 
like interaction strength which exhibits the most strong 
changes as a function of the SP in the itinerant many- 
body system. 

In conclusion, we have determined the real-space e — h 
scattering amplitude and its contribution to the quasi- 
particle self-energy within a local approximation by us- 
ing, as a background, a variational momentum-space 
approach^ We have shown that for low-energy quasi- 
particle excitations this amplitude can be expressed in 
terms of an instantaneous and zero-range spin-dependent 
pseudopotential. The strength V^^, of such a pseudopo- 
tential is defined by the volume integral of the statically 
screened Coulomb interaction with a weighting factor 
concerned with the mutually uncorrelated motion of the 
e — h pair carrying spin 5^ = 0, ±1. The analysis car- 
ried out for the strength shows that the obtained results, 
due to the transparent dependencies on those physical 
variables which are based on Hedin's polarization con- 
cept, could be useful in attempts to go beyond the GW 
method in an ab initio study on itinerant ferromagnets. 
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